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ABSTRACT 



This thesis develops the mathematical expressions 
to evaluate component weight and space for all signifi- 
cant components within the containment vessel (secondary 
shielding) of a Gas Oooled Reactor (GCR) closed cycle 
gas turbine plant for ship propulsion. The equations 
are solved for a range of operating conditions using 
helium as the working fluid. Information from these 
equations was used to obtain the following: 

1. Optimum total weight of components. 

2. Curves of component weight. 

3. Curves of component dimensions. 

The resulting data support the following conclu- 
sions : 



1. Optimum weight occurs at a pressure 
ratio that is lower than the value that gives best 
efficiency. 



2. A heavy weight penalty must be paid 
for small increases in efficiency. 

3. The regenerator is the controlling com- 
ponent in a cycle of this type. As heat transfer re- 
quirements go up (regenerator effectiveness approaches 
100 %), regenerator weight rises sharply and limits 
appreciable gains in efficiency. 

4. Rotating machinery length will have a 
greater influence on plant design than the weight of 
this component. 

The curves should give the designer a rough 
idea of machinery component sizes and weights, and 



ii 



are intended for use in the development and design of 
an efficient, nuclear gas turbine power plant. 



Thesis 
Title i 



Supervisor; S. C. Powell 

Associate Professor of 
Marine Engineering 



iii 



ACKNOWLEDGMENT 



The author wishes to express appreciation for 
the assistance and direction given by Professor S. C. 
Powell, the supervisor of this thesis. 



iv 



TABLE OF CONTENTS 



Page 

Abstract ii 

Acknowledgment iv 

Table of Cont ents . v 

List of Figures vii 

Nomenclature viii 

I. Introduction 1 

II. Procedure 9 

A. General 9 

B„ Heat Exchanger Equipment ................... 14 

C » Reactor 21 

D. Turbo Machinery 26 

E. Piping and Accumulator 30 

F. Evaluation of Constants .................... 35 

III. Results 37 

A. Cycle Conditions, Table III-A 37 

B. Graphs , Table III-B 38 

IV. Discussion of Results 50 

A. General 50 

B. Total Weight. .,.oo.oo..oo.o.oe.oeoo.«ooooo.. 50 

C. Regenerator 50 

D. Precooler 54 

E e Intercooler 55 

F . Reactor o«©. 00000000 . 00 . .. 0000000 . . 0000.0000 55 

G, Rotating Machinery 56 

V. Conclusions 58 

VI . Recommendations 61 

VII. Appendix 63 

A. Heat Exchanger Equations Development 64 

I . General ..••.o.. ©.•©«•©•©•, .©©•.o...... 64 

II . Regenerator, .•...©..©oo©..©...©........ 64 

III. Coolers ............................... 70 

B. Turbo -Machinery Equations Development ...... 73 

I . General .............................. 73 

II. Co mp r essors ....«*. «©o*ooo.oo.©«o*..o.. 73 

III . Turbines 78 

C. Evaluation of Constants 81 

I . General 81 

II. Regenerator..... 82 

III. Coolers 34 

IV. Reactor 87 

V. Turbo Machinery 89 



TABLE OF CONTENTS 



(Continued) 



gage 



D. 



E 



v 



Sample Calculations oo..o.o.©o©oooo.o®ooo©©»« 

1 O 0 en e I* al ooooooooooooooooooooooooo.oo.oo 

II. Sample Calculations 

Literature Citations .00.00000.00.00. ««©.*«. 

I. General Publications 

II. Atomic Energy Commission Publications.. 

III. Special References .............o. 000000 



92 

92 

93 
102 
102 
102 
103 



vi 



LIST OF FIGURES 



Figure 


Title 


Page 


II-l 


Cycle Flow Diagram 


11 


II-2 


Temperature-Entropy Diagram 


11 


II- 3 


Cooler Temperature-Length Diagram 


17 


III-l 


Total Weight and Efficiency 


39 


III-2 


Regenerator Weight 


40 


III-3 


Precooler and Intercooler Weight 


41 


1 1 1-4 


Reactor Weight 


42 


III-5 


Rotating Machinery Weight 


4J 


III-6 


Regenerator Dimensions 


44 


1 1 1-7 


Precooler Dimensions 


4-5 


III-8 


Intercooler Dimensions 


46 


III-9 


Reactor Dimensions 


47 


III-IO 


Rotating Machinery Length and 
LPT Exhaust Diameter 


48 


III-ll 


Flow Rate 


49 


A-l 


Regenerator Temperature-Length Diagram 


64 


B-l 


Stage Velocity Diagram 


77 



vii 
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diameter 
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4(cross section area) 
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heat rate (BTU/HR) 

constant , function of tube arrangement 

gas constant (ft-Lb/Lb°R) 
radius (ft) 

pressure ratio (>l) 

Reynolds number 

entropy (BTU/Lb°F) 

temperature (°R = °F + 460) and °F 

heat transfer coefficient (BTU/ft 2 -HR-°F) 
circumferential velocity on pitch diameter 

K 

exponent = 

velocity (ft/sec) 
volume (ft 3 ) 

K-l 

exponent = — 



work (BTU/Lb) 
weight (tons) 

flow rate (Lb/HR) 

pressure loss factor 




P 

6 

"x 



stress parameter = — — (ft) 

f (° 

angle of relative velocity 

blade stress (p si) 
efficiency 

regenerator effectiveness 
angular velocity (rad/sec) 
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T 


ratio of compressor isentropic exit temperature to 
its inlet temperature 


P 


density (Lb/ft 3 ) or (tons/ft 3 ) 


u 


viscosity (Lb/HR-ft) 


A T 

m 


log mean temperature difference 


Subscripts 


A 


accumulator 


a 


axial 




area 


c 


compressor 

cold 
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diameter 


6 


gas 


HPC 


high pressure compressor 


HPT 


high pressure turbine 


h 


hot 


i 


inside 


ic 


intercooler 


LPC 


low pressure compressor 


LPT 


low pressure turbine 


l 


length 


m 


mean 


N 


number of tubes 


0 


outside 


P 


piping 


pc 


precooler 



X 



R reactor 

RG regenerator 

s surface 

stage 

t turbine 

V reactor vessel 

w water 

we ight 

x cross section 



State Points 



1. precooler outlet 
LPC inlet 

2. LPC exit 
intercooler inlet 

2s. isentropic LPC exit 

3= intercooler outlet 

HPC inlet 

4. HPC exit 

cold side regenerator inlet 

4s. isentropic HPC exit 

5. cold side regenerator outlet 
reactor inlet 

6. reactor outlet 
HPT inlet 

7 . HPT exit 
LPT inlet 

7s. isentropic HPT exit 

8. LPT exit 

hot side regenerator inlet 

8s. isentropic LPT exit 

9. hot side regenerator outlet 
precooler inlet 
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I, INTRODUCTION 



The thermodynamic potential of the closed cycle gas 
turbine power plant has been well recognized for many 
years. Early attempts to effectively utilize this poten- 
tial, however, were unsuccessful because of poor compon- 
ent efficiencies and the lack of a suitable high tempera- 
ture heat source. Technical developments of high tempera 
ture gas turbines over the past years have enabled the 
designer to develop components with the necessary thermal 
efficiency, and the advent of the nuclear reactor has 
created the type of heat source necessary to drive the 
system. 

In the design of a ship propulsion plant, certain 
parameters are not subject to the designer’s control or 
choice. The maximum temperature in the cycle (Tg) may 
be restricted because of reactor fuel cladding creep 
strength limitations or the restriction imposed by the 
high pressure turbine (HPT) blading. The minimum tempera 
ture (T^) of the working fluid is likely to be set from 
the sea water inlet temperature and reasonable size and 
weight coolers. A fixed reactor inlet temperature (T^) 
is necessary to prevent cyclic thermal stresses from 
damaging the reactor structural material and moderator. 

In a gas cooled reactor, the inlet temperature is put at 
a relatively high value to anneal out radiation damage 
in the moderator. Compressor and turbine efficiencies 



- 1 - 



are dependent on the state of technical development of 
these machines. In view of the above, Tg, T^, T^, and 

rotating machinery efficiencies will most likely be 
fixed for any cycle because of mechanical design or 
metallurgical restraints and size of the plant. 

The designer then has at his disposal other para- 
meters which he must fix so as to give the most economic, 
efficient or the least weight and/or space propulsion 
plant . 

The closed cycle nuclear gas turbine plant con- 
sidered in this work consists of two principal sections: 
the reactor, which is located within the primary shield, 
and the propulsion equipment which is located within 
the secondary shielding or containment vessel along 
with the reactor. The main objective of this thesis was 
to develop a weight and space optimizing procedure for 
all significant components located within the contain- 
ment vessel. The machinery components considered as 
significant with respect to the variation in their weights 
and sizes with changes in the cycle parameters are: the 
regenerator, coolers, rotating machinery, reactor, piping, 
and accumulator. The developed expressions for component 
weights and dimensions in terms of thermodynamic proper- 
ties, gas constants and flow rate enables the designer 
to compare the relative size and weight of each component 
for different working fluids. 
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The second thesis objective was to make a weight 
and size calculation of these components using the 
proposed helium-cooled, graphite moderated, Maritime 
Gas-Cooled Reactor (MGCR) power plant as a reference. 
Solution of the developed equations for a range of 
operating conditions provided data for determining 
minimum weight and space for the regenerative, inter- 
cooled gas turbine cycle. 

No gas turbine nuclear power plant has been or 
is presently being constructed for ship propulsion. 
However, with the technological advances being made 
in nuclear reactors, high temperature materials and 
rotating machinery, it is to be expected that a ship 
power plant of this type will someday be used. The 
shortage of component weight and space data for such 
a system will present a major problem to the ship pre- 
liminary designer who must develop an efficient nuclear 
gas turbine plant. The developed weight and space equa- 
tions along with the system study of a particular plant 
undertaken in this thesis should assist the naval archi- 
tect and marine engineer in developing the most economi- 
cal, efficient or the least weight and/or space gas tur- 
bine nuclear power plant. 

For any ship propulsion plant, it will be neces- 
sary for the ship designer to produce machinery compon- 
ents which give the best performance for the least 
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weight, cost, space or some other criteria. The subma- 
rine designer for instance would want to optimize his 
combination of machinery components in order to achieve 
the best performance for the least weight and space. 

The merchant ship designer on the other hand would 
direct his work toward producing a plant that gives 
the best performance for the least cost. Optimization 
of a propulsion plant can be performed for different 
reasons. Knowing the mission, size and requirements 
of the ship, the designer is able to optimize the machin- 
ery on a basis that gives the best power plant commen- 
surate with the ship’s purpose. 

In an optimization procedure, the real problem 
is to reduce the number of variables so that they can 
be conveniently handled. The ideal situation is to 
have the weight or cost expression in terms of one para- 
meter or at most two independent parameters. Computa- 
tion work and time are greatly reduced with each less 
variable. Also desirable is an expression in closed 
analytic form which will greatly simplify calculations. 
However, the functions to be differentiated will usually 
not be found in this form, and the designer must then 
resort to a graphical, numerical or some other method 
of approximate differentiation. It is important that 
careful study be given all parameters in a weight or 
cost function to determine whether or not they are 
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independent. If two parameters are dependent, the method 
of evaluating the function will be different and gener- 
ally more difficult than if all variables are independent. 

The first step in eliminating the number of vari- 
ables is to assign possible values to those quantities 
that are beyond the control of the designer. Secondly 
the concept of ^balance 11 is used to reduce the number 
of variables. The latter step can be applied to compon- 
ents which bear a specific relationship to each other. 

For example, the task of finding the surface area of 
two different components for a given power plant might 
be reduced to a single problem of selecting a balanced 
combination of these areas. This method would be pos- 
sible if the effects of each area on cost or plant 
weight were evidenced by their influences on the same 
working fluid property. The problem of selecting a 
balanced combination of areas is done separate from the 
optimization of the plant as a whole. All other compon- 
ents of the plant must be sensitive only to this fluid 
property and not to the various combinations of these 
two areas which establish that fluid property. Once a 
relationship between balanced components is formulated, 
the units comprising the balanced group are treated as 
a single component with a consequent elimination of 
variables. 

In this work, a number of variables were eliminated 
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from consideration of step one. The ^balance 1 * concept 
was used to form a relationship between the compressors 
and the high pressure turbine. This step eliminated 
the high pressure turbine and low pressure turbine 
pressure ratios as variables. Thus once a reasonable 
value for the cycle pressure loss is specified, the 
remaining variables are flow rate (w) , top pressure 
level (P), overall compressor ratio (r), regenerator 
effectiveness (ti ) and temperatures in the cycle. 

Both cycle efficiency and net work per pound of fluid 
flow are independent of top pressure as shown in Sec- 
tion II. Pressure is then of importance only in the 
way it effects weight and space of components. All 
remaining temperatures of the cycle are dependent on r. 
For a specific power output, flow rate is a function of 
the pressure ratio. With establishment of a fixed reac- 
tor inlet temperature as previously discussed, regenera- 
tor effectiveness is also a function of pressure ratio. 
Because the reactor inlet temperature is set, the cycle 
pressure ratio must always be less than the value that 
would make the turbine exhaust temperature equal to the 
reactor inlet temperature. For equal reactor inlet and 
turbine exhaust temperatures, a ti of 100 per cent is 

A 

required. A pressure ratio giving a turbine exhaust 
temperature lower than reactor inlet temperature would 
require a T| greater than 100 per cent and hence presents 

X 
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an impossible cycle. 

In view of the above, there are two independent 
variables for a plant of this type. All weight and size 
equations were developed in terms of these two inde- 
pendent parameters: pressure and pressure ratio. Sub- 

sequent solution of the equations for a narrow range 
of pressure ratios indicated extreme variations in 
total weight and size of components. 

The regenerator weight and size was found to 
dominate that of all other components by several orders 
of magnitude as regenerator effectiveness approaches 
100 per cent. In order to achieve a higher efficiency 
without having to accept the large increase in weight 
contributed by the regenerator, consideration was given 
to the idea of using an additional piece of heat exchang- 
er equipment between the regenerator and reactor inlet. 

A small quantity of high temperature fluid could be re- 
moved from the reactor to further heat the bulk of the 
fluid from the regenerator. The extraction of this 
fluid from the reactor might occur after it has passed 
up the outside of the reflector in cooling the thermal 
shield and vessel wall and before it flows along the 
surface of the fuel rod assemblies. Another possible 
extraction point for the small quantity of high tempera- 
ture fluid might be in the LPT after partial expansion 

( I \) 

or at the HPT exhaust. This bled fluid would then be 
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recirculated to the additional heat exchanger to further 
raise the temperature of the bulk of the fluid coming 
from the regenerator. A considerable portion of the 
regenerator duty could then be taken over by the addi- 
tional heat exchanger. The resulting lower regenerator 
effectiveness would lower regenerator weight. A study 
of such a proposal is needed to indicate if the reduc- 
tion in regenerator weight, caused by a lower regenera- 
tor effectiveness, is more than the weight increase 
caused by the additional heat exchanger (with associated 
equipment) and the greater weights of other components 
(that might be caused by the slightly larger flow rate 
required in order to keep the same power output). 
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II. PROCEDURE 



A, General 

The basic procedure followed was; 

1. Develop mathematical expressions for component 
weights and sizes in terms of state point properties 
and system constants. 

2, Determine the required system constants. 

5. Solve the equations for a range of operating 
conditions „ 

Steps 1 and 2 are carried out in this section 
with the detailed procedures of the development per- 
formed in Appendices A, B and C. Step 3 is found in 
Appendix D. 

The plant components considered, consist of the 
regenerator, precooler, intercooler, reactor, accumula- 
tor, piping, and rotating machinery comprised of the 
low pressure compressor (LPC), high pressure compressor 
(HPC), high pressure turbine (HPT) and low pressure 
turbine (LPT). The HPT drives the HPC and LPC, and the 
LPT provides shaft power. A low pressure turbine by- 
pass heat sink was not considered as a plant component. 
Preliminary studies [13] indicate that direct bypass- 
ing of the low pressure turbine without an additional 

■^Numbers in brackets refer to references listed in 
Appendix E. 
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heat exchanger is feasible. With this system thermal 
cycling conditions appear to be within the capabilities 
of the different component materials. A system requir- 
ing additional piping, but still no bypass heat sink, 
would use a split-flow scheme. A portion of the by- 
passed fluid would be piped directly to the precooler 
and the remainder would proceed through the regenerator 
as usual. The reactor inlet temperature could be kept 
more closely regulated with this latter method. Figure 
II-l is a flow diagram of the system indicating state 
points and some of the symbols used in the development 
of the equations. Figure II-2 is the corresponding 
temperature-entropy diagram for the cycle. 

The mathematical development that follows is based 
on the perfect gas relationships and on turbulent flow 
of the working fluid in ducts and piping. In addition, 
the velocity of fluid flow in the system is considered 
to be of the order of magnitude where static pressure 
and temperature are approximately equal to stagnation 
conditions. All gas properties are evaluated at the 
mean fluid temperature in a particular component and 
are considered constant over the range of temperatures 
encountered in that component 

For minimum compressor work in a two-stage com- 
pression with intercooling to the low temperature of 
the cycle (T^, the pressure ratio of each compressor 
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~ Transfer 




Figure II/-1 

Flow Diagram for helium-cooled, graphite moderated 
Nuclear Gas Turbine Power Plant 



Figure IIZ-2 

Temperature-Entropy Diagram 
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can be readily shown to equal the square root of the over 
all pressure ratio: 



P 



2 



P 



1 





( 2 - 1 ) 



The above condition also implies negligible pressure 
drop of the fluid through the intercooler. 

Total compressor work can be written for equal 
compressor efficiencies as: 

2c T 

w = (t-5 - 1) (2-2) 

'c 

Total turbine work for equal efficiencies is: 



w t * We (1 - 



( 2 - 5 ) 



Y being the pressure loss factor in the cycle defined as: 



Y = 



^6/^Ss _It ^ M V _ ( V J> ^ ^ h ^T ( 2-4) 

T 4s/ T 1 ‘ V p s P 4 P 2 " V P 5 P 4 P 8 P 9 V C ^ 



t is the ratio of compressor isentropic exit temperature 
to its inlet temperature and is the ratio of turbine 
inlet temperature to its isentropic exit temperature. 

The factors in the last term are the fluid pressure drops 
contributed respectively by the reactor, cold side of 
regenerator, hot side of regenerator, precooler and inter- 
cooler. The pressure drop being a function of working 
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fluid, pressure, and component arrangement and geometry 
is usually estimated by the designer for each component. 
The total cycle pressure drop and hence Y, for the cycle 
or a component, can be determined. 

The net work for the cycle is simply W. - W^. 

Using (2-2) and (2-3) the net work ( w ne ^) P er pound 
of fluid flow can be written as: 



^net _ ^6 /, 

C T, “ n T, K± 

pi 1 



- 1- ^ 
c 



(2=5) 



as : 



Heat into the cycle from the reactor can be written 



Q 



rn 
C T 



P 1 



C p (T 6 - V 

“Vh 



( 2 - 6 ) 



The thermal efficiency of the cycle is then: 



W 

ri = net 



Q. 



m 



(2-7) 



Combining (2-5), (2-6) and (2-7) 



l t T- 



^ " Yt^ 



2 . 

T|, 



(t 0> - 1) 



T r 



( 2 - 8 ) 



T n 



For reasons discussed in Section I, T^, T 1 

t| and tj are not generally subject to the designer’s 

b C 
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control. With a reasonable value assigned to Y, the net 
work and efficiency become functions of cycle pressure 
ratio only. By setting the derivative with respect to 
t equal to zero in (2=5), the optimum value of cycle 
pressure ratio (r) is found which maximizes cycle net 
worko Similarly the optimum value of r can be found 
from (2=8) which gives maximum cycle efficiency,, 

Bo Heat Exchanger Equipment 

All heat exchanger equipment considered is assumed 
to behave as if of the shell and tube counterflow con- 
figuration. In the regenerator the high pressure, cold 
gas flows inside the tube, with the hotter gas flowing 
longitudinally around the tube interspace. The increased 
strength of the tubes with internal pressure would neces- 
sitate this arrangement. Water flows through the tubes 
and gas around the tubes in both coolers. This arrange- 
ment would facilitate the required periodic cleaning of 
the water side of the coolers. The precooler and inter- 
cooler require large cooling water flows with low water 
velocities and low pump heads. Design factors for the 
coolers follow marine steam condenser practice. 

1. Regenerator 

For a given gas, tube size and tube arrangement, 
equations (A-30) and (A-31) can be written respectively 

^Equations carrying letter designations are from the 
corresponding appendix. 
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as ; 



Nd 2 

1 



= K* tj w 



Nd " 1 -t^ 



5 .5 

i ( T c + r “ T h C d> 



(2-9) 



1.4 



Nd i L - K ML * (lirT 



^8 <V r ‘ W* 



( 2 - 10 ) 



where C a = \ (V 

q V o 



s 4.8 



and the K’s are constants 



of proportionality. Since T c and were defined as 



T 



* T 4 



(2-11) 



and = 



T 8 + ^9 



( 2 - 12 ) 



they may be put in terms of pressure ratio and known 
cycle temperatures. From the definition of turbine works 



yT 6 - v - vt (T e - w 



t s ~ v 1 “ y 1 



Yt 



)] 



(2-13) 



Using (2=13) and (A-4), (2=11) and (2=12) can be written: 



V 2 "V 



" T 6 [1 " ^t^ 1 ' Yg^ 

2 o - V 



(2=14) 
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(2-15) 



T. 



h 




5^5 



Substituting (2-14) and (2-15) into (2-9) and 
(2-10), cross sectional and heat transfer area variations 
become functions of presure ratio (r), ti , w, P and 
known quantities: For a specified maximum and mini- 

mum temperature of the working fluid (T^ and T^), a 
fixed reactor inlet temperature (T^) and an appropriate 
value for turbine efficiency, the area variations are 
dependent on r, P, and w. ti becomes known once the 
pressure ratio is specified, since reactor inlet tem- 
perature is fixed. 

Regenerator weight may be written as a function 
of the heat transfer area: 



where K-^ and are constants and (Nd^L) is given by 
(2-10) with (2-14) and (2-15) substituted for T and 

V 

V 

2, Coolers 

For a given gas, tube size and tube arrangement, 
equations (A-32) and (A-33) for the precooler become: 



W RG " K l< Nd i L) ♦ K 2 



(2-16) 




T 



,•5 

pc 



p 



r 



(2-17) 
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.1.4 



Nd o L - K NdL w 



where Zi T = - T-^ 



_£ 



m 



T 



pc 



.4 




,8 



T 



pc 



T 9 + T 1 



(2-18) 

(2-19) 

( 2 - 20 ) 



Using equations (2-13) and. (A-4) , can be expressed as; 



T 6 (1 - 2 "x-"t +2 "t ,| x + 5 t 



2ti. t, 

fH) + n t 



Yt 



T 9 = 



x 5 



( 2 - 21 ) 



1 - T) 



X 



Zl T m , the LMTD may be written as; 



Z\ T = a(zh T. ) = a 
^ m A ' 



/iiT i + AT ( 



7 



= a[T - T ] (2-22) 

pc mw 



where z^.T^ is the arithmetic mean temperature difference 



and u a ; * is dependent on the ratio ^ x i \ see Figure II-3 



AT. 




*AT e 



Figure II-3 

Cooler Temperature-Length Diagram 
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For the precooler and intercooler, a mean water tempera- 
ture (T mw ) of 80°F is assumed, based on 70°F sea water 
inlet temperature and a 20°F water temperature rise 
through each cooler. From [2]^ values of l, a u range from 

1.00 to 0,710 as ^ T i goes from 1 to 10. 

Zi T 

e 

With T w ^ and T 2 specif ied , M a # can be approximated by; 



a = 1.227 - .0006 T 9 (2-25) 

for the range of gas temperatures likely to be encoun- 
tered at the precooler inlet. 

With the above simplifications, (2-19) > (2-20) and 
(2-22) may be expressed in terms of pressure ratio, r\ 

and known temperatures and efficiencies by using (2-21) 
and (2-23). The resulting relations for (2-19), (2=20) 
and (2-22) may then be substituted into (2=17) and (2-18) 
giving rather long equations for the cross sectional and 
heat transfer areas in terms of flow rate, pressure, 
pressure ratio, r) and known quantities. 

For a given working fluid, tube size and tube ar- 
rangement, (A-34) and (A-35) for the intercooler become; 




Nd L 
o 



K NdL w 




1.4 





P 



.8 



(2-24) 



(2-25) 
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where 



(2-26) 



^ T g = T 2 



t 3 = t 2 



T. = 

1C 



T ? : 



T 2 + T, 



(2 



T 2 represents the temperature out of the LPC and can 
be expressed from compressor work as: 



t 2 - T 



-i (t-5 



1 ) + 1 



Combining (2-28) with (2-26) and (2-27): 



AT s = T x 



T-, 






Tl, 



(T^ - 1) 



T ic " t 



1 f -r°5 

*c 

The LMTD corresponds to that of the precooler: 



-1 (t° - 1) + 2 
*1 



AT = a(T. - T ) 
^ m ic nwr 



(S 






(2 



(2 



where' a' ; can be approximated by: 

a = 1.814 - .0014 T 2 (2 

for reasonable values of gas temperature at the inter- 
cooler inlet. Substituting (2-28) into (2-32): 



a = 1.814 - .0014 T, 



— (t* 5 - 1) + 1 



(^ 



-27) 

-28) 

-29) 

-30) 

-3D 

-32) 

-33) 
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Combining equations (2-30), (2-31) and (2-33): 



/\T m = jl.81-.0014 T 1 [^(T‘ 5 -l) + l])^[^(T* 5 -l) + 2]-540j (2-34) 

where T ffiw was taken as 80 + 460 = 54-0°R. 



By placing (2-29) , (2-30) and (2-34) into both (2-24) 
and (2-25), the two intercooler area variations are 
expressed in terms of flow rate, pressure, pressure 
ratio and known temperatures and compressor efficiency. 

As for the regenerator, the precooler and inter- 
cooler weights may be expressed as functions of the 
heat transfer area. For the precooler the basic weight 
equation is: 



MT 

W = K, w/ 

P° 1 zVT m 



\1.4- 



T 



.4 

“pc 



,8 



+ K. 



(2-35) 



where T ,AT and T „ can be expressed in terms of P, 
g 7 in pc 

r, ti x and known quantities from equations (2-19) through 



(2-23). 



The weight equation for the intercooler: 



W. 

1C 



K x w 





( 2 - 36 ) 



where a T , /\ T and T. can be expressed in terms of P, 
g ’ m ic 

r and known values from equations (2-29)* (2-30) and (2-34) 
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C . Reactor 



The reactor weight and size equations are developed 
by considering the reactor as another piece of heat ex- 
changer apparatus. The proposed MGCR helium cooled, 
graphite moderated reactor with a heterogeneous assembly 
of fuel rods [11] , [12] , [14] is used as the reference 
reactor. 

For a constant heat generation flux along the 
coolant channel and radially over the reactor, the gas 
temperature rises with a constant slope versus coolant 
channel length. The coolant wall temperature will also 
rise with the same slope from (T^ + AT) to (T^ + aT) 

where AT is the coolant wall temperature minus the 
temperature of the coolant. The maximum wall and cool- 
ant temperature will occur at the reactor outlet with 
the outlet coolant temperature in each channel being 
equal for uniform heat generation in the core. A con- 
stant heat generation flux is substantially true for a 
relatively small reactor core using good reflector mate- 
rial [1], [3], [15] and [18]. 

The heat transfer in the reactor can be represented 

by: 

Q » H A s AT (2-37) 

Q = w C p (T 6 - T 5 ) (2-37A) 

with A T = maximum wall temperature minus the outlet 



- 21 - 



temperature of the coolant. H is given by (A-13): 



C 1 k ’ 2 (w C p )>8 

N*8 d 1 * 8 (Pr)* 4 



(2-38) 



D and N represent the outside diameter of the fuel rod 
and number of fuel rods respectively. The pressure 
drop through the reactor is given by (A-17). The term 



, /T - T 

1 ( w 

2 \ T 



in (A-17) can not be neglected, but is 



equal to a constant for uniform heat generation and 
fixed inlet and outlet coolant temperatures. (A-17) 
then becomes: 



In ^ = -C 2 ~ [4f (2-39) 

5 

where for low mach numbers, static pressures and tempera- 
tures can be interchanged with stagnation properties. 

The pressure drop factor (Y^) for the reactor is: 



Y 



R 




(2-40) 



Following the procedure used in obtaining equation 
(A-22) in Appendix A, the pressure drop term for the 
reactor can be written: 



y R . C 3 M > i) 



(2-41) 



being a constant determined by the fixed inlet and 



outlet temperatures and the maximum wall temperature of 
the reactor. Using (A-24) and (A-26) the Fanning fric- 
tion factor and mach number squared are represented 



respectively by; 



.2 



f = .0438 



DNC ]i 
C 

w p 



(2-42) 



M 



■Hf fc 



(wC )‘ 
P 



T mR /k-l 1 X2 



RK 



K PND 2 



(2-43) 



mR 



T + T 

65 = constant. Combining all constant terms: 



M 2 = C, 



fe) 2 

KR 



K-l 

K 



PND 2 



(2-44) 



Substitute (2-42) and (2-44) into (2-41): 



. 1. 8 



y R “ C M3 l K 



K-l | ^ (^) L_ (‘"'p) 



y,C ' 



R 



P 2 n-*- * Sj/*’ • ® 



(2-45) 



Combine (2-37) and (2-38) and write: 



C. k’ 2 (wC ) c8 
_i E (tiNDL) 



- N .8 ]) 1.8 ( .p r > ) .4 



(2-46) 



as defined with equations (2-37) is a fixed value 
and so (2-46) can be written: 



G = Q P° 8 ._CPr 

1 k* 2 (wC )’ 8 N c2 L 
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J'- 



(2-47) 



Multiplying equation (2-45) by (2-47) and solving the 
resulting expression for ND 2 : 



ND 2 




(2-57A) can be substituted for Q to give: 



(2-48) 



ND 2 





(Pr)° 3 
R’ 3 



(t 6 -t 5 ) >3 

P 



(2-49) 



The above equation represents the cross sectional area 
variation of the reactor. For fixed inlet and outlet 
temperatures and a given fuel element and coolant channel 
configuration, it is dependent on the flow rate and top 
pressure level in the cycle. 

The corresponding heat transfer area (NDL) equa- 
tion is obtained by solving (2-47) for L, dividing by 
D, making the appropriate substitution for (ND 2 ) * and 
multiplying this result with (2-49); 



NDL 



'NDL 




1.2 



r) 



.4 



..8 



(Pr) 




w (p) (T 6 -T 5 j' 



(£-50) 



In (2-50), L is the coolant channel length or the core 
height, H. Using the bare core approximation, the rela- 
tionship between height and radius for a critical cylinder 
of minimum volume is given by reference [4] as: 



H = 



* ( 2 )' ? 
2.405 



R 



1.85R 



(2-51) 
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D is fixed for a given reactor and N, the number of fuel 
rods (coolant channels) is a function of reactor cross 
sectional area: N = C N R 2 . Placing this and (2-51) into 

( 2 - 50 )^ reactor radius cubed becomes: 



R 5 = C T 






78< T 6-V ' ( 2 -52) 



W R = P R % H- H 



Reactor weight can be expressed as: 

(2-53) 

where is the overall reactor density determined from 
a reference or similar reactor. Using (2-51) and (2-52) 
reactor weight becomes: 

\ 1 « 2 / ^ X . A - /-| X «2 



W R ■ C R 



m (ly ar o. 



^ 8 (V T 5^' 4 (2-54) 



For a given coolant and fuel element configuration it is: 



W R “ K R p. 8 ^ T 6 " T 5^ 



(2-55) 



The R and H used above are for the core and do not 
represent the actual size of the reactor. The overall 
reactor radius and height, exclusive of shielding, can 
be taken as some constant multiple of the core dimen- 
sions for small changes in core R and H. From the MGCR 
reference reactor, the reactor vessel radius (Ry) and 
height (Hy) can be represented as: 

R v = 1.65 00 (2-56) 



- 25 - 



Hy = 3.51(H) 



(2-37) 



D, Turbo Machinery 

1. Compressors 

In formulating length and weight equations for 
the turbines and compressors, it was assumed that the 
length of each component varies directly with the num- 
ber of stages and that weight is proportional to the 
product of annulus area at the low pressure end and 
the length. Therefore by picking reasonable values 
for: efficiency, P, £jj , 3 1 , P 2 and £ , equation 

(B-19) in Appendix B can be written as: 

N LPC = K T i ( T * 5 -D (2-58) 

Since T^, the minimum temperature of the cycle is fixed, 
the LPC length is : 

LlpC " H < T ' 5 - (2-59) 

The annulus area given by (B-12) becomes: 



A LPC 




wr 

P 



(2-60) 



w is the flow rate, r the pressure ratio and P the top 
pressure in the cycle. 

The weight equation, after applying the condition 
stated above to (2-59) and (2-60), is: 
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( t *^ - 1 ) 



(2-61) 



W. 



LPC 




wr 

P 



Using (B-ll) and (B-15), the tip diameter at the LPC 
low pressure end becomes: 



D 



LPC " 





(2-62) 



The four corresponding equations for the HPC are 
indicated below: 



L HPC 

a hpc 

W HPC 



(t- 5 _ i) 



K 



w r 



• 5 



jr E Zll ( r »5 
K w P U 



1 ) 



(2-63) 

(2-64) 

(2-65) 



D HPC ~ K d 





( 2 - 66 ) 



The constants for each of the above equations are 
not necessarily applicable to both the HPC and the LPC. 
Each machine has a particular constant associated with 

it o 



2. Turbines 

As for the compressors, reasonable values of effi- 

ciences, P, £a , p, , and ^ can be assumed or ob- 
U 1 ^ d 

tained from the turbine manuf acturer . These factors will 
change very little over a wide range of operating conditions 
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(pressure, pressure ratio, flow rate, etc.) and so the 
size and weight equations can be expressed similar to 
those for the compressors. 

For the HPT, (B-32) gives the number of stages 

as : 

Njjpp = K (t * 5 - 1) (2-67) 



Turbine length, annulus area, weight and tip diameter 
equations become respectively upon using the basic 
relations derived in Section III of Appendix B: 

- 1) (2-68) 



L HPT ' ^ T< 



A HPT ' K a 



w [1 

P 



TJ t C(f 5 - 1)] 

[1 _ c(t* 5 - l)] u 



2 T, 



where C = 



T 6V. 



(2-69) 



W HPT ' K w 



w (t* 5 - 1)[1 - T|.C(t #5 - 1)] 

p i 

[1 - C(t #!? -1)] u 



(2-70) 



D HPT ~ K d 



I) 



.5 



[1 



T\ t C( T 



.5 



1)] 



.5 



[1 - C( x 



.5 



l )] 1 



(2-71) 



For the LPT, the equations are: 



L LPT = [1 “^t C ^^“ 1)][Y T(l ~ C[T,3 ~ i:i) ~ 1] 

Yt [1 - C(t‘ 5 - 1)] 



(2-72) 
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( 2 - 75 ) 



= K 



wr 



LPT a P 



[1 - r\Al 



Yt 



-)] 



wr [l-tl t (l-Y^][l-n t C(T' 5 -l)][YT(l-C(T- 5 -l])-l] 
W L“ K w P 



Yt[1-C(t' 5 -1)] 



(2 



D LPT ~ K d 



Or )' 5 ^ - V 1 - 



• 5 



(2 



The constants associated with each equation will 
he different for the HPT and LPT unless considerable 
similarity exists between the two machines. 
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E. PIPING AND ACCUMULATOR 



1 . General 

Both piping and the accumulator contribute signifi- 
cantly toward the total weight of the components within 
the containment vessel. In addition, the space occupied 
by the accumulator will effect the size of the contain- 
ment vessel and hence the amount of secondary shielding. 

Weight equations for piping and the accumulator 
are developed in this section. The constants associated 
with these equations were not determined due to lack of 
suitable data from a reference design. The equations 
however indicate the trends to be expected for different 
operating conditions. 

2 . Pip ing 

Weight of piping is considered proportional to the 
internal surface area plus a constant. 



W = K, uDL + K 0 (2-76) 

pi d 



For negligible heat transfer, the pressure drop term for 
the piping is expressed as; 



y 



p 



K-l 

2 



M 2 



L 

D 



(2-77) 



From (A-24-) , 
flow is; 



the Fanning friction factor for turbulent 




(2-78) 



- 30 - 



Mach number squared is given by (A-26) as: 



M 



= C 



(wC p )' 

R 



T_ (K-l) 2 1_ 

T3 2 K s D 4 



(2-79) 



Substitute (2-78) and (2-79) into (2-77) and solve 
resulting expression for D. 



D = K ( ^e!. 041 ( wC )-375 fljS) 

p .41 p V R ) 
p 



• 21 fcl J 62 



\ K 



(2-80) 



For a given working fluid, (2-80) becomes: 



D = K 



w 



• 375 



p .41 

P 



(TL) 



.21 



(2-81) 



Since temperature in equation (2-81) is raised to only 
the 0.21 power, temperature does not effect the diameter 
(weight) significantly. As a preliminary conservative 
weight estimate, the maximum temperature in the cycle is 
used with (2-81). For greater accuracy, the piping length 
could be divided into two halves, a hot and a cold section, 
Equation (2-81) can now be written as: 



w « 375 on 
D - K P (L) 



(2-82) 



Combining (2-76) and (2-82): 



W = K 
P P 



w 



.375 



1.21 



,41 



+ K, 



(2-83) 



Considering the length of piping within the con- 
tainment vessel as a fixed quantity for similar ships, 

(2-83) becomes: 
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(2-84) 



W 

P 



K-, 



w 



• 375 



.41 



K, 



The constants in (2-84) were not determined since 
suitable information from a reference design concerning 
length and weight of piping was not available. 

3. Accumulator 

The accumulator weight is expressed as a function 
of the total weight of the working fluid in the cycle 
and the allowable stress in the wall material. Consider- 
ing a long cylindrical container, the surface area, upon 
neglecting the end plates, is: 



A = 2mRL (2-85) 

The volume of the working fluid is given by: 

V = tcR 2 L (2-86) 

V = m — (2-87) 

P 

where m represents the total mass of the working fluid 
within the cycle. Weight of the accumulator is simply 
the surface area times the product of wall thickness 
and material density, 

W A = 2mRLt (2-88) 

Using simple hoop stress analogy, wall thickness becomes: 

t = ( P ~ F a ) R (2-89) 
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where P a is the external (atmospheric) pressure and P 
is the pressure of the fluid within the accumulator. 

For rapid changes of load, the pressure of the fluid 
within the accumulator must be greater than the top 
pressure of the cycle. Equation (2-87) indicates a 
high pressure is also desirable from the standpoint 
of reducing the size of the container. To facilitate 
rapid changes in plant load [10] , the accumulator is 
arranged in two parts : a high and low pressure sec- 

tion. When reduction in loads occur, gas is bled from 
the HPC discharge to a low pressure accumulator or 
receiver. As the low pressure flask fills, it is 
emptied by a piston type compressor which compresses 
the gas to a pressure considerably above the top press- 
ure in the cycle. The small piston compressor discharges 
to the HP accumulator. As loads again increase, the HP 
accumulator discharges into the high pressure line 
between the HPC and the regenerator. Equating (2-86) 
and (2-87) and solving for the radius squared (R 2 ): 

D 2 _ mRT 

K ” tcLP 

Combining (2-88), (2-89) and (2-90), the weight of 
accumulator becomes: 

w 2 mRT ( P ” P a ) 

W = 

A Op 



(2-90) 



(2-91) 






Since the temperature of the gas to be put in the ac- 
cumulator will be of the order of 200 - 300°F (the HPC 
discharge gas temperature), the accumulator material 
need only be a carbon steel. Therefore for a given 
gas and a moderate storage temperature, equation 
(2-91) reduces to: 

W A = K m (P ~ P a) (2-92) 

6 P 

Because the pressure within the accumulator will be 

several orders of magnitude larger than atmospheric 
P-P 

pressure, ^ and (2-92) becomes: 

P 

W A = K | (2-93) 

Maintaining a constant hoop stress for a fixed cylinder 
radius means that the pressure to wall thickness ratio 
must be a constant. With these restrictions, accumulator 
weight is written: 

W A K A m (2-94) 
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F. Evaluation of Const ants 



The mathematical development of the weight and 
space equations necessitated the introduction of a num- 
ber of unknown constants. Some constants depend upon 
materials used, some are proportionality constants and 
others are constants of design. 

To solve the equations, it is necessary to evalu- 
ate each constant. Accurate evaluation of the constants 
requires a thorough investigation of similar plants, 
with note taken of the specifications, weight state- 
ments, heat balances and other information. 

The equations developed under Section II and 
Appendices A and B are general in nature and not re- 
stricted to any particular working fluid or type of 
ship. Only in the evaluation and application of the 
constants does the problem become one of a specific 
type or design. 

Since no shipboard nuclear closed cycle gas tur- 
bine plant exists, available data for evaluating the 
constants was limited. The author relied mainly on 
information from the proposed MGCR propulsion plant 
presented in references [11], [12], [14], [15] and, 
where applicable on information from existing open 
cycle gas turbine plants. The heat balance, specifica- 
tions, and weight estimates for the helium-cooled, 
graphite moderated MGCR plant allowed satisfactory 
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evaluation of the constants for all components except 
piping and the accumulator. 

The evaluation of these constants and sample 
calculations are given respectively in Appendices C 
and D. 
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III. RESULTS 



A. Cycle Conditions 

Original results were obtained in the form of data 
from hand calculations of the developed equations with 
reference to the helium cooled, graphite moderated MGCR 
power plant. Calculations were made by varying the 
compressor pressure ratio (r) for two different reactor 
inlet temperatures (T^): = 700°F and T^ = 79l°F. 

Table III-A lists the fixed values used for both sets 
of calculations. 



Table III-A 

Fixed Cycle Conditions 



t 

LPT (net power) output =2l,500HP Pressure d!rop=— =.86 

Working Fluid =Helium Coolant Water Inlet 

Temperature(T w ^) =70°F 

Top p ressure (P) =1000 psia Coolant Water Outlet 

TemperatureCT^) =90°F 

Maximum Cycle Temperature Turbine Efficiency 

(T 6 ) =1300°F (rj t ) =.90 

Minimum Cycle Temperature Compressor Efficiency 

(T x ) = T 5 = 100°F (n c ) =.86 
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B. Graphs 

The data is presented in graphical form in accordance 
with Table III-B. 



Table III-B 
Graphs 

Plot Figure 

Total Weight and Efficiency vs Pressure Ratio III-l 

Regenerator Weight vs Pressure Ratio III-2 

Precooler and Intercooler Weight vs Pressure Ratio III-3 

Reactor Weight vs Pressure Ratio III-4 

Rotating Machinery Weight vs Pressure Ratio III-5 

Regenerator Dimensions vs Pressure Ratio III-6 

Precooler Dimensions vs Pressure Ratio III-7 

Intercooler Dimensions vs Pressure Ratio III-8 

Reactor Dimensions vs Pressure Ratio III-9 

Rotating Machinery Length and Low Pressure 

Turbine Exhaust Diameter vs Pressure Ratio III-10 

Flow Rate (lb/HP-HR) vs Pressure Ratio III-ll 
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PRESSURE RATIO 



IV, DISCUSSION OP RESULTS 



A. General 

The results obtained in this thesis are based on 
the Maritime Gas-Cooled Reactor (MGCR) plant and are ap- 
plicable only to a plant of this type. Accuracy of re- 
sults are dependent on the accuracy with which system 
constants were evaluated. 

Hand calculations were performed to solve the 
developed equations for a range of operating conditions. 
Results are presented graphically in Section III. 

For all components, the sharp increase in weight 
or cross sectional area indicated by the curves at low 
pressure ratios is a result of the large flow rates re- 
quired at these ratios. The pressure ratios which give 
a mathematical regenerator effectiveness of 100 per cent 
are approximately 3»1 and 3»8 respectively for the high 
and low reactor inlet temperatures. A higher pressure 
ratio than 3*1 for the cycle with T^ = 791°F would pre- 
sent an impossible plant since a regenerator effective- 
ness greater than 100 < /^o is required. A similar state- 
ment applies to the cycle with the lower reactor inlet 
temperature. 

B. Total Weight 

Total weight and efficiency are plotted against 
pressure ratio for two different reactor inlet 
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temperatures in Figure III-l. Minimum weight occurs at 
a pressure ratio of about 2.5 and gives an efficiency- 
equal to 30.5 °Z> for T^ = 791 °F. For a reactor inlet 
temperature of 700°F a slightly higher minimum weight 
occurs at a pressure ratio equal to 3.15 with efficiency 
again equal to 30»5 ^o. For pressure ratios greater 
than that for optimum weight, the weight rises sharply 
to infinity due to the regenerator weight. This is be- 
cause of the large increase in regenerator heat transfer 
area required as turbine exhaust temperature approaches 
reactor inlet temperature. The MGCR proposed plant is 
designed for a pressure ratio of 2.6, a value slightly 
above that for optimum weight, but giving an efficiency 
of 32 %. 

Figure III-l indicates that a greater gain in ef- 
ficiency is possible with small changes in pressure ratio 
above that ratio which gives optimum weight. Although 
lowering the reactor inlet temperature gives a flatter 
optimum weight, overall efficiency is significantly re- 
duced since more heat must be put into the cycle for the 
same reactor outlet temperature. In addition, the opti- 
mum weight for T^ = 700°F is larger than that for a 
higher inlet temperature because of the increased reac- 
tor weight. This greater weight is caused by the increase 
in temperature rise associated with a low inlet temperature 
as indicated by equation (2-55)* The lower inlet 
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temperature would, also tend to add to the complexity of 
reactor construction. 

C. Regenerator 

Regenerator weight reaches an optimum at a lower 
pressure ratio than that for total weight, and rapidly 
approaches infinity as tj goes to 100 c / 7 b. Prom Figures 
III-l and III-2, it is noted that the regenerator is the 
controlling element in total weight at high pressure 
ratios. The designer desires a high pressure ratio in 
order to achieve more efficiency, hut without the large 
increase in weight that must he accepted as indicated 
hy these figures. 

A 100 % regenerator effectiveness can never he 
realized, and in reality the designer must he willing 
to pay heavily in terms of weight and space for a regen- 
erator effectiveness as low as 90 °7o . When one con- 
siders the weight and efficiency curve for = 791°F» 
it is worthy to note that efficiency can he increased 
5 points by going from a pressure ratio of 2.5 to 3»0* 
However, for this increase in efficiency, total weight 
is more than doubled. Undoubtedly this sizeable gain 
in efficiency could not be realized for such a weight 
penalty. 

A possible solution in obtaining part of this in- 
crease in efficiency with a less drastic increase in 
weight might be to use an additional piece of heat 
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exchanger equipment between the regenerator and reactor 
inlet. A small quantity of high temperature working 
fluid could be taken from the reactor. This extraction 
might occur after the fluid has passed up the outside 
of the reflector in cooling the thermal shield and ves- 
sel wall and before it flows along the surface of the 
fuel rod assemblies. Another extraction point for the 
small quantity of high temperature fluid might be in 
the LPT after partial expansion or at the exhaust of 
the HPT. This extracted or ff bled W fluid would then be 
recirculated to the additional heat exchanger to further 
raise the temperature of the bulk of the fluid coming 
from the regenerator. A considerable portion of the 
original regenerator duty would then be taken over by 
the additional heat exchanger (pre-heater). The result- 
ing lower regenerator effectiveness would lower regenera- 
tor weight. The heat transmission in this proposed pre- 
heater would be very high since both sides of the addi- 
tional heat exchanger will be near the top pressure in 
the cycle. With the high rate of heat transmission pos- 
sible under these conditions, it may be possible to keep 
the size of the pre-heater quite small. Although this 
small quantity of bled fluid would again enter the 
reactor at a higher temperature than the bulk of the 
gas, the thermal inertia of the reactor and the large 
quantity of fluid should prevent any harmful cyclic 
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thermal stresses. In order to produce the same power, 
a larger flow rate of fluid is required with the addi- 
tional heat exchanger. In addition, a small pump 
(compressor) would be required to raise the pressure 
of the *\led' J gas to reactor inlet pressure. A study 
of such a proposal would indicate if the reduction in 
regenerator weight, caused by a lower regenerator effec- 
tiveness, is more than the weight increase caused by the 
additional heat exchanger (with associated equipment) 
and the greater weight of other components associated 
with the slightly larger flow rate. Plant complexity 
with this system may be such that would not warrant its 
use. 

Figure III-6 indicates the variation of regenera- 
tor length from zero to infinity as pressure ratio in- 
creases and more heat transfer area is required. The 
curves indicating length in this figure are for each 
unit of the 3 unit regenerator component . 

D. Precooler 

Figure III-3 shows for a pressure ratio of 2 and 
larger that precooler weight has a slight negative 
slope. This is to be expected since there is less and 
less heat transfer area required as the turbine exhaust 
temperature is lowered with increasing pressure ratio. 
For the lower reactor inlet temperature, the precooler 
weight increases. This effect is due to the increased 
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heat transfer area required of the precooler as regenera- 
tor effectiveness decreases. 

Precooler length rises with a decreasing pressure 
ratio as shown in Figure III-7. This would he neces- 
sary since there is little temperature drop across the 
turbines and the fluid must be essentially cooled from 
the reactor outlet temperature to the fixed low tem- 
perature of the cycle. A significant increase in 
length is required for T^ = 700°F since the precooler 
must perform more heat transfer as tj decreases. 

X 

E. Intercooler 

The intercooler weight, as shown by Figure III-3, 
increases slightly at higher pressure ratios. As the 
LPC discharge temperature is raised with pressure ratio, 
more heat transfer area is required to reduce this tem- 
perature to the fixed low temperature of the cycle. 

There is no change in intercooler heat transfer area 
(weight) with variation in reactor inlet temperature. 

The intercooler must always cool the gas from the LPC 
discharge temperature to the cycle’s bottom temperature. 

Figure III-8 indicates that intercooler length in- 
creases with pressure ratio. For low ratios, little 
heat transfer area is required to bring the gas back to 
the cycle’s minimum fluid temperature. 

F. Reactor 



Reactor weight increases sharply at low pressure 
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ratios since a bulky heat generating unit is required 
to raise the temperature of the large flow rates. 

Since reactor weight is directly proportional to flow 
rate, weight will decrease with less flow. It is also 
worthy of note that for a drop of 91 °F in reactor in- 
let temperature, a weight increase of approximately 
35 tons results at a pressure ratio of 3»0* This is 
to be expected since reactor weight as given by (2-55) 
is proportional to temperature rise to the 1.4 power. 
Reactor vessel diameter and height curves have the same 
form as the reactor weight curves. The lower inlet tem- 
perature requires larger reactor dimensions since equa- 
tion ( 2 - 52 ) shows that radius cubed is proportional to 
the temperature rise in the reactor to the 1.4 power. 

G. Rotating Machinery 

As indicated by Figure III-5, rotating machinery 
has a minimum weight for a pressure ratio of slightly 
less than 2.0. At low pressure ratios, the large an- 
nulus area required causes an increase in the weight. 

At higher pressure ratios, the weight increase is not 
as sharp since length of turbo machinery has less ef- 
fect on weight than the diameter or annulus area. 

Rotating machinery length being a strong func- 
tion of the number of stages also increases sharply 
with pressure ratio. Few stages are required at low 
ratios and many stages are needed at higher pressure 
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ratios 



Tip diameter for only the low pressure turbine 
exhaust is indicated in Figure III-10 since its dimen- 
sions will he larger than that of the other 3 rotating 
machines . 

Although the rotating machinery weight is sig- 
nificant, the length of the unit is likely to have a 
greater influence on design of the plant. A machine 
of great length would necessitate a larger contain- 
ment vessel with a corresponding increase in shielding 
weight. An in-line arrangement of compressors and tur- 
bines would probably result in lower total rotating 
machinery weight and fewer shaft sealing, bearing and 
support problems than an arrangement where the low pres- 
sure turbine is in a separate casing. The latter ar- 
rangement would permit a shorter overall length by 
placing the LPT beside the HPT and compressor unit. 

This arrangement would have a slight detrimental effect 
on efficiency because of losses in piping between the 
HPT and LPT, but may be necessary if containment space 
is limited. 
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V. CONCLUSIONS 



An analysis of the results supports the following 
conclusions : 

1. Optimum weight for the reference plant occurs at a 
value of pressure ratio that is lower than the value 
that gives best efficiency. 

Efficiency at optimum weight (pressure ratio of 
2.5) is 30.5 °7o. This compares with a maximum possible 
efficiency of 35*5 at a pressure ratio of 3»1 

(corresponding to a regenerator effectiveness of 100 ^q) 
for a reactor inlet temperature of 791 °F» 

For a reactor inlet temperature of 700°F, the ef- 
ficiency at optimum weight (pressure ratio of 3 * 15 ) is 
30.5 %. The maximum possible efficiency is 33*3 °7o 
occurring at a pressure ratio of approximately 3»8 
(corresponding to a regenerator effectiveness of 100 ^o). 

2. With a plant of this type, the designer must be will- 
ing to accept a heavy weight penalty for slight effici- 
ency gains as pressure ratio is increased above the 
optimum weight value to the ratio that corresponds to a 
regenerator effectiveness of 100 ^o. Table V-A indi- 
cates the changes in efficiency and total weight for 
small increments of pressure ratio above that for 
optimum weight. 
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Table V-A 



Changes in Efficiency and Total Weight vs Pressure Ratio 



Pressure ratio increment Zsn (°7o) 

Reactor inlet temperature 
2.5 - 2.6 1.0 



2.6 - 2.7 1.0 

2.7 - 2.8 0.8 

2.8 - 2.9 0.8 

2.9 - 3.0 0.9 



W(tons) 

791°F 

5 

7 

15 

55 

290 



Reactor 

3.15 - 5-25 
5.25 - 5.55 
5.55 - 5.45 
5.45 - 5.55 
5.55 - 5.65 
3.65 - 5.75 



inlet temperature = 700°F 

.5 4 
.4 5 
.4 10 
.5 16 
.4 40 
.5 165 



3. A lower reactor inlet temperature gives a flatter but 
higher optimum weight. A significant rise in reactor 
weight results from a lower reactor inlet temperature. 

4. The regenerator is the controlling component in 

this cycle. As heat transfer requirements go up (regen- 
erator effectiveness approaches 100 c /'o), regenerator 
weight rises sharply and limits the gain in efficiency 
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that is otherwise possible with a larger pressure ratio. 

5. The precooler and intercooler weights and dimensions 
experience a minimum over a wide range of possible 
operating pressure ratios and do not strongly influence 
the total weight. 

6. Length restrictions imposed by available space are 
much more likely to influence rotating machinery design 
than the weight of this component. 
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VI. RECOMMENDATIONS 



The following recommendations are made: 

1. Since most of the developed equations show that weight 
is inversely proportional to pressure (to a power), the 
tendency would be to increase pressure level to reduce 
weight and plant size. However, maximum pressure limi- 
tations would be set from the high temperature , pressure 
effects on materials. For a more complete analysis an 
investigation should be made to determine the effect of 
pressure on the reactor, piping and regenerator weight. 

2. Equations should be refined to give greater accuracy. 
Accuracy can be improved as plants of this type are de- 
signed and built. A survey of these plants for compon- 
ent weight and size data will contribute to a more pre- 
cise evaluation of system constants. Such a survey 
would provide information for proper evaluation of the 
piping and accumulator constants. 

3. Shielding accounts for a major portion of total 
plant weight, with the secondary shielding weight being 
largely a function of size and arrangement of components 
within the containment vessel. A thorough study of size 
and arrangement of components is necessary to ensure a 
minimum shielding weight. Reference [93 gives a suit- 
able method for the preliminary evaluation of shielding 
weight . 
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4. The regenerator weight and size dominates that of all 
other components by several orders of magnitude as regen- 
erator effectiveness approaches 100 *7*0 . An investigation 
as to the feasibility of recirculating a small quantity 
of fluid from the reactor after it has cooled the thermal 
shield or from the turbine (after partial expansion) through 
an additional heat exchanger (pre-heater), located between 
the regenerator and reactor inlet, is necessary. Such a 
study would indicate if the reduction in regenerator weight, 
caused by a lower regenerator effectiveness, is more than 
the weight increase caused by the additional heat exchanger 
and the greater weight of other components associated with 
the slightly larger flow rates. 
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APPENDIX A 



HEAT EXCHANGER EQUATIONS DEVELOPMENT 

I. General 

Development of equations for the heat exchanger 
cross sectional and heat transfer area variation are 
based on reference [17]. 

Pressure drop of the gas due to entry and exit 
effect in the exchangers are considered negligible. 

II. Regenerator 




2 - distance along regenerator 
Figure A~1 

Regenerator Temperature-Length Diagram 

The rate of heat transfer from the hot fluid to 
the cold fluid for equal flow rates and heat capacities 
on each side may be written in the form: 



Q = 


wC p 


(T 5 - V 


(A-l) 


Q = 


wC 

P 


< T S ' V 


(A-2) 


Q = 


UA . 

s 


AT 

m 


(a-5) 
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where AT m = = xxT g is the log mean temperature dif- 

ference across the counter flow regenerator shown in Figure 

A-l . 



Regenerator effectiveness is defined as: 
-l„ = *5 - ** - T 8 - 



T 8 - T 4 



T 8 - T 4 



Combining equations (A-l), (A-2), (A- 3) and (A-4) 



V ■ at: 



w C 



m 



(hj 

' ‘x 



(A-4) 



(A-5) 



The overall coefficient of heat transfer is written as 



UA 



H~J — + l n(d o/ d i ) + 1 

n h sh k^lrcL 



H A 
c sc 



(A-6) 



For a gas-to-gas heat exchanger, the thermal resis- 
tance of the tube wall may be neglected in comparison with 
the film resistances and (A-6) may be written: 



1-1 

U “ H, A . + 

h sh 



H, 



i / H 

1 / T C 

= H 'c ( + H h 



d i 



(A-7) 



where A = A 

s sc 



c c \ n o ' 

, the heat transfer area inside the tube. 

Reference [ 5 ] expresses the heat transfer coefficient 

of a turbulent fluid as: 

H = .023 -p (Re)’ 8 (Pr)’ 4 (A-8) 

e 

The hydraulic diameter outside of the tubes (hot 
side) becomes 



d , = d q 

eh o 4 - 



where q for triangular tube spacing is: 

ZL&l /b_V _ x 
* 1<V 



q = 



(A-9) 



(A-9A) 



The equivalent outside flow diameter per tube can be 
written as: 

d x = d o (q) ’ 5 (A-10) 
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For the cold side, Reynolds number can be written: 



Re) 



and for the hot side: 



Re), 



4w C 

,td i ^C°p 



4w C 



(A-ll) 



nd Nji, C 
o 'n p 

Substituting (A-ll) into (A-8) and noting that 



(A-12) 



d. = d:- for the cold side: 

Q 1 



8 



IF 



,028 k c* 

N° 8 d. 1<,S (Pr) 

1 c 



(A~13) 



(A-8) becomes for the hot side: 

.2 



H, 



.028 (wC p ) 



.8 



N 



,8 , 1.8 



(A-14) 



q(PrJ 



o 1 '" / h 

Combining equations (A-7)» (A- 13 ) and (A-14) 



1 

U 



N“ 8 d, lo8 (Pr) A 
i c 



1 + 






(A-15) 



where = 



.028 (WC ) * 8 k ° 2 
p c 

^ c ° 2 ( Pr ) ° ^ 

\ 1 for the temperature range 



(Pr) 



.4 



considered. The overall heat transfer coefficient then 
becomes : 



U = .028(WC p )° 8 k c - / 1 ) (A-16) 

N .8 di 1.8( Pr ) c «4 \i + ° 8 / 

a i 

The pressure drop of a gas flowing through a tube 
with combined friction and heat transfer can be expressed 
for low maah numbers [?] as: 
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dP KM 2 / 

--o 



4f d 



1+ I 



T - T 
w o' 



o “ ' " ' ~o 

where f is the Fanning friction factor defined as 

shear stress 

/V/2g 



d(sj (A-17) 



In the regenerator and coolers, the 



o / (p — T \ 

term if _2£ o) is negligible compared to 1 and 

2 \ T / 



is 



dropped. With the low mach number likely to be encountered 
in the ducting and heat exchanger equipment , static pres- 
sures and temperatures can be substituted for stagnation 
pressures and temperatures with little loss of accuracy. 
With these simplifications, (A-17) can be integrated to 
give for the cold side : 



M 



2 n ^ = -K — § 4 f 



d i 



(A-18) 



and for the hot side: 



^ J? ■ 



M. 



= -K 



4 f v 



V 



(A-19) 



From equation (2-4), the pressure drop factor 
(Y^q) contributed by the regenerator is: 



L RG 



\ P 8 



1 



( A-20) 



Take antilog of (A-18) and (A-19) and combine with 
(A-20) to obtain: 

/ v M w f v^TT 

exP (- | (A-21) 

is a constant for the working fluid used and for 



-RG 



a particular top pressure in the cycle. 

After taking the logarithm of (A-21): 
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y RG " 



+ 



(A-22) 



K-l 



)n 



M 2 f 
c c 

d. 

1 



IV 

V / 



y^G i s a constant since Yj^,, as noted above, is constant, 

The relationship for the Fanning friction factor 
(f) for turbulent flow can be given as: 



f 



,046 

(Re)' 2 



Substituting (A-ll) and (A-12) into (A-23): 



f c = .044 



d i N ^n \ 

wG p ) 



.2 



c 



f h = * 044 




.2 

h 



(A-23) 



(A-24) 



(A-25) 



From the mach number definition and the continuity 
relation, M 2 for the cold and hot side may be written 
respectively as: 



M 2 = 



M L - 



4J X 2 

71 



(wc y 
e n p 



fc 



RK 



RK 



K-l 1 
K P Nd( 



' c 



K-l 



K P h Nqd; 



where T 

c 



T 8 + T _2 



T 5 + T 4 and T, = 

2 h 2 

Substituting (A-24), (A-25), (A-26) and (A-27) 

(A-22) and simplifying, the following expression is 

obtained : 



(A-26) 

(A-27) 



into 
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(A-28) 



y RG " C RG 




L 

R 



l 2 i 

n 1.8 ,4.8 -P* 



1 



T +r 2 
c 





where it was assumed (uC ) (uC ), at the mean 

pc p h 

temperatures of the cold and hot side, and all con- 
stants have been written as Cg^. The pressure term 
in the denominator of (A-28) is the top pressure level 
in the cycle. 

Combining (A-3), (A-5) and (A-16): 



wC 



'x 



P 1-th 



.028 (wC )’ 8 k * 2 

V c 

N° 8 d i 1,8 (Pr) 

c 



A + Q 



fc) 



,8i" Nd i L 



(A-29) 



Divide equation (A-28) by (A-29) and solve the 
resulting expression for Nd^s 




Equation (A-30) represents the cross sectional area 
variation of the regenerator in terms of gas properties, 
tube size and arrangement, regenerator effectiveness, 
flow rate, top pressure level, cycle pressure ratio, and 
mean temperatures on the hot and cold side. represents 

the constant terms. 
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The heat transfer area variation (Nd^L) is obtained 
as a function of the same properties by solving (A-29) for 
L/d^, making the appropriate substitution for (Nd^)*2 
and multiplying this result by (A-JO): 



It is noted that equations (A-30) and (A-31) simplify 
considerably for a given gas, tube size and tube arrange- 
ment. The resulting expressions then become functions of 
flow rate, regenerator effectiveness, top pressure level, 
pressure ratio and mean temperatures on the hot and cold 
side. In the last two equations, the gas properties may 
be evaluated at the mean temperature on the hot or cold 
side since Pr and the ratio of specific heats can be 
considered constant over a large range of temperatures 
for most gases. Also the change in thermal conductivity 
to the 0.2 power would vary only slightly since k under- 
goes little change with limited temperature variation. 

Ill . Coolers 

Development of the corresponding cross sectional 
and heat transfer area equations was done using the 
same analogy as for the regenerator. For the coolers 




4 
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1 

I 



I 









the heat transfer coefficient was based on the outside 
of the tubes (gas side). With the approximation that 
the tube wall and water film thermal resistances are 
small compared to the gas thermal resistance, the over- 
all heat transfer coefficient was based solely on the 
gas heat transfer coefficient. 

For the precooler, the cross sectional area vari- 
ation is then given by: 



Ndl = C 



Nd 



w /^ T g 



AT 



m 









iPri 






(A-32) 



The heat transfer area variation is: 



Nd o L =°IBL 



w 






Similarly for the intercooler, the following results 



Nd o " C Nd w ( AT 



/K-lY~*^ (Pr)^ rl 5 



1C 

E / \ K 



(A-34) 



Nd o L ■ C NdL W U\T 



f T icV 



a. 2 



m 



R 



¥y 



(Pr) 



.64 r o 



da }* 2 

-t ^ C«5) 



where R, K, k, and Pr are the familiar gas properties. In 
the above expressions, the gas properties are evaluated at 
the mean gas temperature in the respective cooler. 



T 



T Q + ^ 



pc 



and T 



T^ + _ T 2 + T 1 



1C 
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P and r are again the top pressure level and pressure 
ratio respectively of the cycle. For the precooler: 



A T 



S 




T 



1 



and A T m is log mean temperature difference across 
the precooler. 

For the intercooler: 



AT = T~ - T, = T 0 
g 2 3 2 



T 



1 



andAT„ is the LMTD across the intercooler, 
m 
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APPENDIX B 



TURBO -MACHINERY EQUATIONS DEVELOPMENT 

I. General 

The assumptions involved in the development of the 
rotating machinery equations include; 

1. Uniform axial velocity over the length of the 
component . 

2. All machines of free vortex design. 

3. Equal work done per turbine or compressor stage. 

4. Static pressure and temperature equal to stag- 
nation pressure and temperature. 

5. Steady flow. 



Development of the rotating machinery equations 
are in part based on reference [6]. 



II. Compressors 

The root stress in a blade caused by centrifugal 
force can be expressed in accordance with references 
[6] and [8] as; 



6 = f 



f. -dL 

g 2 

& o 



(£) 



(B~l) 



where f is the blade taper factor and Z 7 the density of 
blade material. Defining a stress parameter (0) as 
equal to * (B-l) is written as; 



0 = 



2g o 




(B-2) 
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Since the annular flow area is: 
A = n&/ 



(B-3) 



(B-2) can be written as: 



P - 



Ji 



2 t S 0 2 S 0 

u 






(B-4) 



where the relation = $J2 was use< ^ obtain the 

term on the right in (B-4). From continuity and the 
perfect gas relationship: 

wRT 



AC a - 



(B-5) 



60 SV 



Using (B-4), (B-5) and the relation RPM 
machinery rotational speed can be expressed: 

, K. K - 5 

RPM = 30 





.5 


[2 s 0 ’ n 


71 w'RT 


t 


H 7 i) J 



(B-6) 

From (B-6) it is seen that RPM may be limited because of 
stress considerations. For a fluid with low sonic velocity, 
RPM will be limited by mach number considerations. For 
RPM designed with mach number considerations, equation 
(B-6) becomes: 



RPM = 



U 



l 



.5 






.5 



(B-7) 



( K S 0 RT ) 8 5 

after substituting (B-4) for All remaining turbo 

machinery relationships will be based on stress considera- 
tions. For helium as a working fluid, compressor and 
turbine design would be influenced by stress considerations 
rather than mach number limitations. 



(B-3) can be written as 



. i 

A = Ttd jp 



(B-8) 
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and 



~ ( nT7d ) 



(B-9) 

Using (B-4-) and (B-5), the pitch diameter is represented 
by: 



d = 



2 wRT 



• 5 r f 

2 s 0 F d 



•.25 



(B-10) 



For the LPC inlet, (B-10) becomes: 

5 / 



d LPC - K 



wT-^r 



'a/U' 



• 5 H 



-.25 



(B-ll) 



where P and r are the top pressure level and the pressure 
ratio respectively of the cycle, and K represents the 
constant terms. 

Substituting (B-ll) into (B-8), the inlet annulus 



area for the LPC is : 



A LPC ' K 



w’T-j^r 



O) 



-.5 



a) 



• 5 



7 v d 

1 



(B-12) 



For most compressors, P, ^ , and can be assigned 

fixed values over a range of operating conditions. Then 
pitch diameter and annulus area are dependent on the 
volumetric rate of flow and cycle pressure ratio. 

Using the same procedure, the corresponding equa- 
tions for the HPC are: 



d HPC ' K| 



w'T 1 ' * ^ 

~T 

wT- 



(r) 



25 



A HPC ~ K p 



^) ,5 | c 



'a/ W 

1 \ 



5 / 0 \ -.25 



a/ u, 
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(ft 
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(B-13) 



(B-14) 



The tip diameter at the inlet end of the LPC can be 
expressed as : 

^ d 



\ - 



*({) 



LPC 



(B-15) 
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with a similar expression for the HPC. 

Applying Newton's law in the form of angular momentum, 
along with the assumptions given in section I of this ap- 
pendix, the work done per compressor stage is: 



W 



“ “§ ( tan ^2 " tan ^ 



(B-16) 



where ^ an< ^ are the angles as shown in Figure B-l, 



For the LPC, the work is given by: 



W 



LPC 



J C T, 
E_1 

tv. 



(t° - 1) 



(B-17) 



Dividing (B-17) by (B-16) the number of stages is: 
JC t> T 1 (t ’ 5 " 



N. 



LPC 



u 2 ® 

c '' sc g~ “u ( tan " tan ^ 



ti t 



(B-18) 



where ti is the stage efficiency. (B-18) is in the 
s c 

form whereby mach number considerations might dictate 
the number of stages since 



. ke So t 1 . 

U 2 KRg Q U* KRg^Lr 

By substituting for U 2 from equation (B— 4), 

be put in the form: 

2JCT, (t° 5 - 1) (/ /d), 

N = 

Vsc* -§ ( tan ^2 " tan ^ 



(B-18) may 



(B-19) 



In (B-19), stress considerations now influence the number 
of stages. 

For the HPC an expression similar to (B-19) is 
applicable since the inlet temperature to both compressors 
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is equal, and only the ^ ratio and the angles an< ^ 
may vary slightly from that of the LPC. 





Figure B-l 

Stage Velocity Diagram 
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III. Turbines 



Before developing the corresponding equations for 
the turbines, it is necessary to determine the pressure 
ratio across each turbine as a function of the overall 
cycle pressure ratio. Since the HpT drives only the 
two compressors, the work of the HPT must equal the 
total compressor work: 



W HPT " ^tS^e " T 7 s^ " l 1 x 



HPT 



(B-20) 



From equation (2-2): 



*c ■ 



2 C T-. c 

P 1 (t* 5 - 1) 



(B-21) 



Equating (B-20) and (B-21) and solving for x 



HPT* 



C HPT 



1 - 2T 1 



(x‘ 



1) 1 - C(t° - 1) 



(B-22) 



where C is a constant since T^, T^, and T) c will be 
fixed for a given cycle. 

From equation (2-4) : 

^ t HPT^ t LPT^ = T t = Yt (B-23) 

Substituting (B-22) into (B-23): 



t lpt = Yt ^ 1 “ 



5 _ i M (B-24) 

At the HPT exit (LPT inlet) the temperature can be 
expressed as a function of T^ and pressure ratio. Using 
the definition of turbine work and (B-22): 

T ? = T 6 [1 - Ti t C(T‘ 5 - 1)] (B-25) 

Following the procedure used for the compressors 
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and the preceding work in this section, the pitch dia- 



meter and annulus area at the HPT exit can be written: 



d HPT ” K| 



2T, 



/WT 6 [l-ri t C(T- 5 -l)] 
P[1 -C(t’ 5 -1)] u 



• 5 






'a/U, 






(B-26) 



where C = 



t 6 ¥< 



HPT 



wT,[1-ti C(t )5 -1)] c (£\.'? 

- K — ^ K O 



P[1-C(t' 5 -1)] U 



'a/U 



(B-27) 



For the low pressure turbine: 



/wT, 



d LEP~ Kl 



A LPT" K 



Cl— T1+. ( 1— y^) ]r ) ,5 /l \ 



w ' r 6 [1 -’ 1 t (l -Y?~ )]r 



'a/U/ 



'a/U/ 



• 5 



1 f*3 



O) 



^)- 2 5 
• 5 



-.5A1 



(B-28) 



(B-29) 



The tip diameter at the exhaust end of the LPT is 
given by: 



IL = 1 + t 



'8 



LPT 



(B-30) 



For the HPT, the number of stages is: 

JOpT, (l- 

b V T fnpp / 



N 



HPT 



Vst 



(B-31) 



U 



— jj - (tan ^ 2 - tan ^) 



With (B-22) substituted for (B-31) becomes 



N- 



HPT 



2JCpT 6 c( T ^-i)n t n st U\ 

o ^a (tan - tan P-,) d v 

P U ^ A 



(B-32) 



where (B-4) was used to substitute for U 2 in (B-31); C 
is again the term defined with (B-22)„ 



For the LPT, the number of stages become: 
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! 



N 



IPT 



J V? 



1 - 



t lpt 



Vst 



II 2 

§ (tan - tan P 1 ) 



(B-35) 



Using (B-4), (B-24), and (B-25 )t 



2JC ^Cl-^C(T-?-l)][YTa.CtT- 7 -l])-l]l1 t%t (d/ fi (B _ 54) 



'1 



N. 



IPT" 



P “§(tan ^ - tan P 1 )Yt[1-C(t° ;? -1)] 



The turbo-machinery equations developed in this and 
the preceding section appear to be quite cumbersome. How- 
ever, it should be noted that many of the factors (effi- 
0 C 

ciencies , g-, P, P 2 and p^) will be fixed by the tur- 
bine and compressor design, and will be applicable to a 
wide range of operating conditions, '^hen the only vari- 
ables left to the marine engineer’s discretion are Tg, 

T-p t or r, Y and volumetric flow rate. 



- 80 - 



APPENDIX C 



EVALUATION OF CONSTANTS 



I . General 

The MGCR reference plant has a reactor rating of 
49.7 MW (thermal) with a flow rate of 7^.6 lb/sec. 

The coolant top and low temperatures in the cycle are 
respectively: 

T 6 = 1300°P and ^ = 100°F. 

Reactor inlet temperature is T^ = 79l°£’. The ratio 
of turbine pressure ratio to compressor pressure ratio 
is 0.86, giving a pressure drop of 140 psi through the 
cycle. Compressor pressure ratio is 2*6 with the top 
pressure of the fluid at 1000 psia. Turbine and com- 
pressor efficiencies are respectively: 

t|. = 0.90 and n = .86. 

X/ c 

With the above conditions, the LPT produces 20,600 HP # 
and the HPT 7^6 auxiliary horsepower in addition to 
driving the compressor unit. In this work, the HPT 
provides only sufficient power to overcome the com- 
pressor work, leaving the LPT with an output of 21,500 
HP for the above conditions. 

The remaining temperatures of the reference cycle 
are given below. The state points correspond to those 
indicated in Figure II-l. 



11 

VM 

11 


100°P 


Tg « 852°F 


T 2 * T 4 ■ 


239°F 


T^ = 301 °F 


T 5 - 


? 91 °F 


T wl - ?0°F 


T 6 - 


1300 °F 


T w2 = 90°F 


To = 


1008°F 





II. Regenerator 

The regenerator is composed of three equivalent 
units arranged in series. Since no heavy internal 
supporting structures are necessary, the regenerator 
can he constructed in a U-shaped fashion to conserve 
space when a high regenerator effectiveness is needed. 
This shape unit would also simplify the thermal expan- 
sion problems. The constants evaluated for the regenera- 
tor dimensions and weight are based on the overall size 
and weight of the regenerator. Unit length and weight 
can then be found by dividing by 3» Each unit uses 
3/8 inch outside diameter tubes, 20BWG, with a 5/8 
inch center-to-center triangular tube spacing. Heat 
transfer and cross section area and effective tube length 
are : 

NdL = 23,340 ft 2 
Nd 2 =7.1 ft 2 
L = 81 ft 
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Total weight of the reference regenerator having an 
effectiveness of 0.90 is: 



WT = 37.I tons 

For the above tube data, q as given by equation (A-9A) 
in Appendix A is equal to 2.06. With the above infor- 
mation, equation (2-9) is written: 

. 5 

Nd ' - k nd w (^) p- < T C + • 0421 W V 5 (c-D 

with T c and T^ defined by (2-11) and (2-12) respectively. 
Using (C-l) and the above information is evaluated 
as : 

K nd = .91 (C-2) 

Upon using (2-10), becomes: 

k NDL - 22 5-5 (C- 5 ) 

Number of tubes is given by: 



N = 



cross sectional area 

? u ♦ 1) 



Nd 2 

l d o U - 1) 



(C-4) 



where the denominator equals 2.35 x 10 ^ for the given 
tube data. 

Effective tube length can be found by using: 



heat transfer area _ NdL 

n N d it N d 

o o 



(C-5) 



Total regenerator weight as given by (2-16) is: 
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W RG = (Ndli) + K 2 

where NdL is equation (2-10). Assuming headers and sup- 
port equipment weigh 3 tons per unit, K 2 is taken ass 

K 2 = 9.0 tons (C-6) 

Using (2-16), (C-3) and the given operating con- 
ditions, is then evaluated ass 

K 1 = 1.21 x 10-3 tons/ft 2 (C-7) 

Cross sectional and heat transfer area and the 
regenerator weight can then he written ass 

» 5 

Nd 2 = .9l(w)fe- V ~ (T c + ,042l(r) 2 T h )* 5 (C-8) 

\ ' T l x / 

/Tj \1.4 , .4 

NdL = 223.5 wM-j L 8 (T c + .042l(r) 2 T h ) (C-9) 

x P * 

W RQ = 1.21 x 10“ 5 (NdL) + 9 (C-10) 

with NdL given by (C-9). 

Ill . Coolers 

Both coolers are similar in arrangement and con- 
struction. Each uses 5/8 inch outside diameter tubes, 

18 BWG, with a 7/8 inch center-to-center triangular tube 
spacing. Weight, effective tube length and cross sec- 
tional and heat transfer area for the precooler used in 
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the reference plant are respectively: 



W = 13.6 tons 
pc 

L = 19.2 ft 
Nd* = 11.6 ft 2 
NdL = 7600 ft 2 

Using equation (2-17), the above information and 
the conditions stated in part I of this appendix, 
is evaluated as 

K nd = 1.48 (C-ll) 

With (2-18) , becomes: 

K tol = 290 (0-12) 

Number of tubes and length are given by expressions simi- 
lar to (C-4) and (C-5) respectively. 

Precooler weight is found by using (2-35) ; 



V - k i w Ut £ 



1.4 



m 



pc) 



.8 



+ K 0 



(2=35) 



with 4 T and T given by (2-19) and. (2-20) respectively, 
g pc 

AT ffl is the log mean temperature difference. 



Assuming a header weight of 5 tons : 



K 2 = 5 tons (C-13) 

Using (2-35), (C-12), (C-13) and the given conditions, 
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I 












K-^ is evaluated ass 



= 1.13 x 10 ^ tons/ft 2 



(C-14) 



This value of is of the same order of magnitude as 
that for a steam condenser [ 9 ] but somewhat smaller. 

Weight, effective tube length and cross sectional 
and heat transfer area for the intercooler of the ref- 
erence design ares 

W. = 13.4 tons 

1C 

L = 14.0 ft 

Nd 2 = 13.4 ft 2 

NdL = 6690 ft 2 

With (2-24) and the given conditions s 

% D - 3*21 (C-15) 

Given conditions applied to (2-25) evaluates K NDL ass 



Equation 



K NDL 

(2-36) 



= 477 

gives intercooler weight ass 




K l w 





(C-16) 



(2-36) 



K 2 for the intercooler is considered larger than that 
for the precooler because of the higher pressure associ- 
ated with the intercooler. 

K 2 = 6.5 tons (0-17) 
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Combining (2-36), (C~16), (C-17) and given conditions: 



K 



= 1.03 x 10-3 tons/ft‘ 



(C-18) 



Summarizing, the equations for the precooler and 
intercooler ares 

Precooler: 

•5 



Nd = 1.48 w 



V 



5 r 

pc p 



NdL = 290 V 4 (i ]' 8 

W = 1.13 x 10 ” 5 (NdL) + 5 
pc 



(c-19) 

(C-20) 

(C-21) 



with NdL given by (C-2 0). 



Intercooler! 



Nd 2 = 3.21 w 



NdL = 477 w 



,T_ 
iAT 



• 5 



_£ 



m ■ 



,T 



1.4 



T 



m 



m • 5 E 

1C P 



ip 

ic — q 

P O O 



W. * 1.03 x 10^3 (NdL) + 6.5 



ic 



with NdL given by (C-23) 



(C-22) 

(C-23) 



(C-24) 



IV. Reacto r 

Core radius cubed for the reference reactor was 
developed in sub— section C of Section II and can be 
written from equation ( 2 - 52 ) as: 
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w 



(C-25) 



R 3 



K 



w 

R p • 8 





1.4 



With the core diameter of the MGCR reactor equal to 6.87 
ft, and for the given operating conditions, Kg is evaluated 
from (C-25) as; 

K r = .0176 (C-26) 



The weight of the reference ractor is approximately 
173 tons. Combining (2-51) and (2-53) « 



W R = 5.81 R 5 (C-27) 

Substituting (C-25) and (C-26) into (C-27) and 
applying the given conditions, reactor density is evalu- 
ated as; 

(° R = .92 tons/(ft 3 of reactor core) (C-28) 

Considering the reactor vessel (exclusive of 
shielding), = .105 tons/ft 5 . Substituting (C-25), 

(C-26) and (C-28) into (C-27), reactor weight becomes; 



w R = .0945 (t 6 - t 5 ) 1,4 ^3 (c-29) 

Core radius cubed is given by; 

R 3 - .0176 (T 6 -T 5 ) le4 (C-50) 

Core height can then be found using (2-51); 



H = 1.85R 
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(2-51) 



Reactor vessel radius (Ry) and height H.y are given 
by (2-56) and (2-57): 

Ry = 1.65 (R) (2-56) 

Hy = 3.51 (H) (2-57) 

V. Turbo Machinery 

In developing the turbo machinery equations, a 
large number of variables are present. Fortunately most 
of these variables are a function of the machine design 
and not a strong function of the cycle operating condi- 
tions. Therefore most similar machine variables will 
have values that are near those indicated below (for the 
reference plant) for a considerable range of operating 
conditions . 

Both machines : 

P 10,000 - 12,000 ft 




Compressor Turbine 

(low pressure end) .11 .20 

efficiency .86 .90 

Weight, machine length and tip diameter of the 
rotating machinery are considered to be of significance 
in determining the overall plant size and weight. There- 
fore the constants associated with each of these three 
developed equations (weight, length and tip diameter) were 
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evaluated for each rotating machine. 

The following data apply to the machinery used in 
the reference design; 





LPC 


HPC 


HPT 


LPT 


L (ft) 


7.95 


10.23 


3.98 


2.84 


W (tons) 


4.8 


3.92 


3.42 


3.33 


D (inches) 


18.1 


14.3 


20.2 


25.0 



The applicable length (L) , weight (W) and tip dia- 
meter (D) equations from sub-section D of Section II were 
used in conjunction with the above given geometry data 
and the operating conditions to evaluate the constants. 



The following 


values were obtained 


as constants: 








LPC 


HPC 


HPT 


LPT 




K -e 


(ft) 


37.0 


47.6 


18.5 


27.2 




K w 


Lbf sec 
in 2 


113.0 


151.0 


156.0 


216.0 




K a 


( sec Lbf V 5 

y Lbm J 


40.7 


41.1 


63.4 


65.3 






For the 


LPC, the 


working equations are; 








L = 


37«0(t° 5 


- 1) 






(C-31) 




W = 


113.0 fF 


(t- 5 - 1) 






(0=32) 




D = 


40.7 (f^ 


)- 5 






(C-33) 




For the 


HPC; 












L = 


4?.6(t° 5 


- 1) 






(C-34) 
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(C-35) 



W = 151.0 ~ r’ 5 (t-5 _ 1 ) 

D-^l.l^fJ * 5 r ’ 25 

For the HPT; 

L = 18.5 (t‘ 5 - 1) 

W = 156 = (T- 5 -i)[i-y(T- ? -i)] 
[1 - C(x° 5 - l)] u 

D = 63 . 4^|)° 5 
2T 1 K 

where ° = and u = icH 

For the LPT: 

27.2 [1-T). C(t° 5 -1)][Yt(1-C[t o5 -1])-1] 

_ L: 

Yt [1-C(t° 5 -1)] 



fl-t| t C(T o5 ~l)] 
[1 _ C(t’ 5 -1)] 



(C-36) 

(C-37) 

(C-38) 



(C-39) 



(C-40) 



W = 216.0 



[l-tl t (l-Y^)Hl-n t C(T° 5 »l)][YT(l-C[T o5 -]])-l] (C 
Yt[1 - C(t° 5 - 1)] 



D 



65.3 




(1 




(C-42) 



“41) 
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APPENDIX D 



SAMPLE CALCULATIONS 



I. General 

In order to indicate the effect of certain para- 
meters as they individually influence plant efficiency, 
weight and space, it is desirable to fix as many other 
variables as possible. As stated in Section I, the 
top and low coolant temperatures, reactor inlet tempera- 
ture, turbo-machinery efficiencies and pressure losses 
for the cycle will generally be fixed for a power plant 
of this type. From consideration of the developed equa- 
tions, the remaining parameters effecting the cycle are 
pressure (P) , flow rate (w) , pressure ratio (r), and 
regenerator effectiveness (t x )« However, as noted in 
Section I, regenerator effectiveness is specified once 
a pressure ratio is selected for this type cycle. 

Solution of all equations was performed for a 
LPT (net power) output of 21,500 HP with the below fixed 
conditions s 



T 6 = 1760°R 
T x = 560°R 

T^ = 1251°R 



V = «90 



c 



86 



P = 1000 psia 



t wl " 70 ° F 



r t 

— = .86 or Y 

r 



o 



94-1 




“92 




















I 






Another set of calculations were made with the 
reactor inlet temperature lowered to 1160°R, with all 
other parameters above held constant. 

Equations were evaluated for pressure ratios rang- 
ing from 1.3 to a ratio that placed the LPT exhaust 
temperature (Tq) at the same level as the reactor in- 
let temperature (T^). This condition would require a 
regenerator effectiveness (ti ) of 100 °7o* A n y pressure 
ratio larger than this value would present an impossible 
cycle since tj would have to be greater than 100 °/o . 

For equations where w and P appear, w is in 
lb/sec and P is in p sia. All temperatures are in 
degrees absolute (°R). 

II. Sample Calculations 

Sample calculations are indicated for all equations 
using a pressure ratio of 2.4 and the fixed conditions 
given in Section I of this appendix. 

A. Flow rate 

Flow rate is found using equation (2-5) • 



net power = w [rj^ 



+ W 1 - & 



? T i c „ 



(t- 5 - 1)] 



90(1760)1 . 25 



'(o 941) (1.42) 



J„2£ssxn2a (L19 . 1)] 



w = 81.1 lb/sec 
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